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Global monochromatic solutions of the scalar wave equation are obtained in flat 
wormholes of dimensions 2+1 and 3+1. The solutions are in the form of infinite series 
involving cylindirical and spherical wave functions and they are elucidated by the 
multiple scattering method. Explicit solutions for some limiting cases are illustrated 
as well. The results presented in this work constitute instances of solutions of the 
■ scalar wave equation in a spacetime admitting closed timelike curves. 
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^ ■ I. INTRODUCTION 

\0 ; 

y—^ \ In this paper global solutions of the scalar wave equation in 2+1 and 3+1 dimensional flat 

\ wormholes are obtained. As a significant consequence of their non-trivial topology, worm- 
^ ■ holes admit closed timelike curves (CTC's). As such they constitute a suitable framework 
O ' for the study of the solutions of the scalar wave equation in a spacetime admitting closed 
^ \ timelike curves. Again due to the topology of a wormhole, no single coordinate chart is 
1/-^ \ sufficient to express the global geometry of the whole wormhole spacetime and it becomes 
O ■ necessary to develop techniques to handle global issues on the one hand and to investigate 
O . the propagation of scalar waves near closed timelike curves. It should be mentioned that 
^\ there are works that study the scalar waves that are valid locally in a certain regions (such 
^1 as may be termed the "throat" ) of the wormhole [1] . 

• Wormholes are widely studied and discussed, especially after the paper of Morris and 

. Thorne, in the context of time travel [2], [3]. The wave equation attracts attention in this 
rS I sense that whether the anomalies of causality violations due to CTC's has corresponding 
' footprints in the solutions of wave equation. 

The metric of the spacetime, as well as the topology, maybe the origin of the closed 
timelike curves, as in the case of Godel's universe [4]. Recently Bachelot has studied the 
properties of wave equation on a class of spacetimes of this type [5]. It may be an interesting 
question, whether the origin of the CTC's (being topology induced or metric induced) effects 
the global properties of the solutions. Nevertheless, as proposed in the concluding remarks, 
the presence of CTC's due to topology does not seem to have a significant effect on the 
solutions of wave equation. 

Cauchy problem of the scalar wave equation in the flat wormhole considered here is 
studied throughout by Friedman and Morris with a variety of other spacetimes admitting 
closed timelike curves [6], [7]. They also proved that there exist a unique solution of Cauchy 
problem for a class of spacetimes, including our case, with initial data given at past null 
inflnity [8]. 

Due to the wormhole structure, the boundary conditions imposed in solving the Helmholtz 
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equation depends on the frequency. Therefore spectral theorem is not apphcable in a straight 
forward manner to express the solution of the wave equation as a superposition of monochro- 
matic wave solutions found in this work. However, in [8], it is proved using limiting absorp- 
tion method that, the superpositions of the monochromatic wave solutions of the problem 
converge to the solution of wave equation. 

The problem can be handled as a Cauchy problem with given initial data at past null 
infinity or alternatively as a scattering problem, i.e. finding scattered waves from wormhole 
handle given incident wave. 

Our approach is similar to that used in scattering from infinite parallel cylinders [9]. 
and ^2 represents outgoing cylindrical (or for 3-1-1 dimensions spherical) waves emerging 
from the first and the second wormhole mouth respectively. In order to be able to apply the 
boundary conditions conveniently which arise from the peculiar topology of the wormhole 
in our case, it is necessary to express in terms of cylindrical (spherical) waves centered at 
second mouth and vice versa. Addition theorems for cylindrical and spherical wave functions 
are employed for this purpose. 

The equations for the scattering coefficients of ^'i and ^'2 that result from the boundary 
conditions in question are in general not amenable to direct algebraic manipulation . The 
multiple scattering method is applied to obtain an infinite series solution. On the other 
hand for some important limiting cases the equations solved explicitly. The solutions by 
these both methods are consistent with one other. 

The outline of the paper is as follows: In section II, the spacetime is described and the 
general formulation of the problem is presented. In section III, 2-1-1 dimensional case is 
studied. The equations are presented, explicit solutions for two limiting cases are obtained, 
and finally the multiple scattering solution is applied. In Section IV the same scheme as 
section III is followed for 3+1 dimensional case. In section V numerical verifications of the 
results obtained in section III are presented. Section VI contains some concluding remarks. 



II. FLAT WORMHOLE 

Given a Riemannian manifold M, a solution F : M x R — > C of the scalar wave equation 



is said to be a monochromatic solution with angular frequency a; e R — {0} if it is of the 
form F{m,t) — ^(m)e*'^* for some ^ : M — > C. Clearly ^ is a solution of the Helmholtz 
equation 

A^ + cj2^ = 0. (1) 

On a general Lorentzian spacetime the concept of monochromatic solution makes sense 
provided the spacetime has an almost product structure that singles out the time direction 
locally. 

A simple example of a wormhole topology is the flat wormhole described in [8]. This 
3+1 dimensional flat wormhole spacetime is constructed as follows: Let a,d,T e M with 
d > 2a > 0. Consider 
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FIG. 1: 2+1 dimensional flat wormhole. P is identified with Q. Arrows indicate the direction of 
the identification. 

where A_|_, A_ are open balls of radius a > and respective centers (0, 0, d/2), (0, 0, —d/2) in 
M^. The boundaries of A_|_ and A_ are designated as 2+ and S_ respectively. The wormhole 
spacetime M of width d, radius a, and lag r is the semiriemannian manifold obtained as the 
quotient space of M by identifying events P, Q on E_,_ x R, E_ x R respectively if P is the 
reflection of Q in the xyt-plane after a translation by r along the t-axis, the semiriemannian 
metric being naturally inherited from the ordinary Minkowski metric on M^. M is clearly a 
flat Lorentzian spacetime. To be precise: 

E+ = {{x, y, z) e R'V^ + y^ + {z- d/2f = a^}, 
E_ = {{x, y, z) e R^\x^ + y^ + {z + d/2f = a^}. 

For {x, y, z) e E+, P and Q are identified where 

P = {x,y,z,t), 

Q = {x,y, -z,t + T). 

In 2+1 dimensions the manifold is defined in the same way except that: 

N= (r2-(A+UA_)), 

A_|_, A_ are open disks of radius a > with respective centers {d/2, 0), {—d/2, 0) in R^ and 
P is the refiection of Q in the yt-plane after a translation by r along the t-axis. 
The geometry for 2+1 dimensions is shown in fig. 1. 

Two wormhole conditions arise from this identification map defining the topology. These 
conditions will function as boundary conditions imposed on the general solution of Helmholtz 
equation in a flat spacetime. 

The two wormhole conditions will be denoted as C-1 and C-2. C-1 is 



F{P) = F{Q), 



FIG. 2: Coordinates used for 2+1 dimensions. 



and C-2 is 

np-VF{P) = -nQ-VF{Q). 

where uq is the unit outward normal to E_ at Q and np is the unit outward normal to S_|_ 
at P. In terms of C-1 and C-2 are: 

np ■ V^{uj,p) = -e'^'^fiQ ■ V^icu, q), 

where p and q are the projections of P and Q on N respectively. 

The solution will be expressed in three components: An everywhere regular part of the 
wave, which may be considered as originating from the sources at past null infinity (or 
alternatively as the incident wave if the problem is considered as a scattering problem) , and 
two outgoing waves originating from each wormhole mouth (or scattered waves from each 
mouth), ^1 and ^2- Obviously ^ = + + ^2 

There are two wormhole conditions that enable one to determine two of \E'o, ^E'l and ^2- 
The problem will be handled like a scattering problem and the scattered waves ^'i and ^'2. 
will be solved given the incident wave ^'o- 

III. 2+1 DIMENSIONS: 

In 2+1 dimensions, solution of Helmholtz equation in cylindrical coordinates yields Bessel 
(or Hankcl) functions. Being everywhere regular, \l/o is expressed in terms of J„(r), while 
^'land ^'2 represent outgoing waves radiating from the wormhole mouths A_ and A_|_, re- 
spectively. Outgoing waves are expressed by Hankel functions of the first kind, Hn\r). 
Referring to fig. 2, has its natural coordinates {r,0) centered at (— (i/2,0), and ^2 has 
its natural coordinates (i?, 0) centered at ((i/2,0). The coordinate variables, 9 and (j) are 
chosen in this way to make use of the mirror symmetry of the geometry of the wormhole 
with respect to y axis. Since ^'o, ^^land ^'2 are valid in exterior domain, they are expressed 
in terms of integer order Bessel (Hankel) functions only. Therefore the expansion of ^^o, 
^land ^2 in terms of Bessel (Hankel) functions are: 
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n=— oo 
oo 

n=— cx) 
oo 

*2= J] CnH^^\uR) ■ e'^t 



Bn and C„ will be found given the coefficients of the incident wave An- The two wormhole 
conditions supply the two equations to determine the unknown coefficients and C„. 
The wormhole conditions C-1 and C-2 are: 



d d 

— ^|jj=a,0=9 = -e^^ — '^\r=a,e - TT < ^ < TT. 

oR or 

To compute ^' at i? = a and r = a it is necessary to write down ^'q, ^'i in {R, 4>) 
coordinates and in {f^, G) coordinates. The addition theorem for cylindrical harmonics 
is used for expressing a cylindrical wave in terms of cylindrical waves of a translated origin 
[10]. It should be noted that, unlike the everywhere regular Bessel functions Jn(a;r), there 
are two different versions of the addition theorems of Hankel functions. For f — d + R, 
addition theorems yield 



E Jkiood)H^^l^{ur)e'^''+''^'^ if r > d 

fe=— oo 

E Jkiujd)Hl^li^{ujR)e-'^^+''^'f' a R>d 
//«(a;r)e-(^-)=<| , (3) 

E Hl^\u;d)Jn+k{ujR)e-'(''+''^^ if R < d 

k=—oo 
oo 

J„(a;r)e^"(^-^) = ^ Mu;d)Jn+k{coR)e'^''+''^t (4) 

fe=— oo 

Wormhole conditions require the expression at r = a and R — a. Since a < d, r < d 
versions of (2) and (3) should be used. 

Accordingly, the wave functions are expressed as a sum of Bessel functions at translated 



origin as 

oo oo 

n=— oo n=—oo 
oo oo 

n=— oo n=— oo 

oo oo 

n=—oo n=— oo 

The expressions for An, and Cn are found using (2), (3) and (4) as: 

oo 

Y ^k-nJk{ujd), (5) 

fe=— oo 
oo 

Bn^ Yl Bk-nHi'\u;d), (6) 

fe=— oo 

oo 

Cn= Y Ck-nH'i^\ujd). (7) 

fc=— oo 

Having obtained the expression of the wave in the coordinates centered at each mouth, 
apphcation of wormhole conditions give necessary equations for the unknown coefficients B^ 
and Cn- 

C-1 leads to 

oo 

Y {^n ■ Jni^a) + Bn ■ H!^\uja) + CnJn{^a))e''^' 

n=—oo 

oo 



and C-2 leads to 



Bn ~ S^^^Cn — — ^^(1) ~{^n ~ C^'^^^n + Cn — e*'^^i?„), (8) 



oo / o a f9 \ 

n,=—oo ^ ^ 

oo / o \ 

Y [AnQ-^U00r)\r=a + S„ — + C n^H^^\ur)\r=a\ e^"', 

1— — rv-i ' 



OO 



— e 

n=— oo 
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_dr 
d_ 

dr 



+ e^--C„ = {A^ + e^-M„ + + e^^^E^). (9) 



Solving (8) and (9) for _B„ and C„, one finds 

S„ = -^+{ua)Cn + e^'^^7^(^«)^n - 7n (^a)^n + e^'^TrT (^a)^, (10) 
C„ = -7+(c^a)S„ + e-^'^^7-(a;a)C'„ - 7rt(^a)^ + e-'^^%{uja)A^, (11) 

where 



7n(^«) 



In (^«) - 2 



or 



For the sake of simplicity the known parts of (10) and (11) will be denoted by En and 
Fn respectively. 



En = -jU^a)An + e"'^%iua)An, (12) 
Fn = —f^{ua)An + e-'^^%{uja)An. (13) 

This pair of equations (10) and (11) are not solvable explicitly; however it is possible to 
solve them for the limiting cases a <^ d and a <til. 



A. Solutions for a <^ d and a <C 1 : 



The difficulty in solving (10) and (11) arises from the convolution sum present in the 
expressions of Bn and Cn- However this term can be evaluated for special forms of Hn\uid), 
namely when it is in complex exponential e~*"" form. When \n\ <^ ud, asymptotically 

Hn\ojd) becomes e"*"""/-^ as a function of n. 7^(a;a) and 7^ (a;a) are almost zero for |n| > 
2a;a, and so are B^ and C„. Thus when a <^ d, the only terms that contribute to 7^(cc;a)5„ 
(7^(a;a)C'n) are those satisfy \n\ < 2u!a <S cud. The a <^ d case is of practical importance 
in physics. In a wormhole universe, this corresponds to the case that the wormhole is 
connecting regions of the universe that are spatially far from each other compared to the 
radius of the wormhole. 

This approximation is not valid for the high frequency limit in general. 

When a ^ 1, 7^(a;a) tends to zero unless n 7^ 0, regardless of d. Accordingly, so are S„ 
and Cn- 

These two cases in which approximate solutions are possible, a <^ d and a -C 1, are 
examined below. 
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1. a<^d 

For large cuo?, asymptotic formula for Hn^ (uod) is 



in-K/2 



(14) 



,^ An^-l .2(4n2-l)(4n2-9) JAv?-l)iAn'^-<d)iAv?-2h) 



l!(8cud) 



2!(8a;(i)^ 



3!(8c^d)3 



2;(a;(i) 



^i(a;d-(7r/4)) 



7ra;(i 



For 'n? <^ the infinite sum inside the brackets can be approximated to 1 : 



This form of Hn\ojd) allows one to evaluate the sum B^- 



fe=— oo m=—oo m=— oo 

= ^(a;d)S(7r/2)e-^"^/' 
where hat denotes the Fourier sum: 

oo 
m=— oo 

Substituting into (10) and (11) 



S„ = z(a;d) [-(7(77/2)7+ + e"^^ Bi7r/2)%]e-'^-/' + E„, 



(15) 
(16) 



The right hand sides of (15) and (16) involves C{7r/2) and B{ti/2) which are unknown 
yet. Multiplying each side by e~*"'^/^ and sum over n's gives a pair of equations for C{v:j2) 
and S(7r/2) : 

(17) 

The numerical results comparing the solutions obtained by these formulae and by the 
multiple scattering method is presented in the appendix. 

4n^ — 1 

To have a better approximation, the second term ^7777^ — 77 in the infinite sum of in (14) 



-B{n/2j 




C{^/2)_ 





'I- z{ujd)e^'^^^-{n) 


z{(jjd)^^{Tl) 


-1 


'E 


>/2)" 


z{ujd)'^^{T{) 


1 - ^(cjrf)e^*^^7-(^) 




F 


>/2). 



can be included: 



l\{%ujd) 



H^:\ud) 



^ iM-(7r/4)) -m^/2Q , i 

nud 1!(82) 
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In this case, the expression for Hn (tod) involves n^e~'"'^/^. This form of H^\u!d) still 
enables one to evaluate i?„ and On, explicitly. This time there arise 6 unknowns and (17) is 
replaced by a 6x6 matrix equation. In this way it is possible to have better approximations 
by taking more terms into account in (14). The number of the linear equations is 4/c — 2 
when the first k term is taken into account in (14). 



2. a < 1 : 

At ua = 0, the Bessel function J„(a;a) is a discrete delta function with respect to variable 
n, and its derivative is zero for all n : 



= \ otherwise ' |:'^'^('^^)l-o = ^ ^■ 



(18) 



Therefore, in the limit a goes to zero, 7^(a;a) and 7„ (cua) become discrete delta functions: 



7+(^a) ^ 7„ (uja) 



where [11], 



^o(^a) 



1 



-Sn, 



To 



l + -(ln(^^)+ 0.5772)' 
n 2 



(19) 



The 7^(c<;a) factors standing in front of each term on the right hand sides of (10) and 
(11) make S„ and C„ delta functions as well. 

Bn = BqSji, 
Cn = CoSn] 

SO that. 



Bq and Co are found by substitution to (10) and (11): 



7±(a;a)S„ ^ {^o5n){BoH^\ujd)) = joBoHi'\u;d)5n, 



'Bo 









1 



- -foHo{u;d)e'^- 7o^o(^c?) 
-foHoiud) 1 - ^oHo{u;d)e-'^\ 



-1 


Eq 







where 



10 



m=— oo 

oo 
m=— oo 

B. Multiple scattering 

An alternative approach is the use of multiple scattering [9], [12]. In the multiple scattering 
method, the scattered waves are decomposed into lower order scattered waves from each 
mouth. Initially each wormhole mouth is considered to be excited by only the incident wave 
and first order scattering coefficients are found by imposing wormhole conditions. Then each 
wormhole is considered to be excited by only the first order scattered wave from the other 
mouth and the second order scattering coefficient are found imposing wormhole conditions. 
f^th q]^(Jqj. scattering coefficients are found by continuing the same procedure. The scattered 
wave from each mouth is the sum of these /c*'^ order scattering coefiicients. 



oo 

k 



k=l 



where 



Yl B'^Hi'\u;r) ■ e^^ and S„ = ^ 



k=l 



Similarly, 



*2 = J]*^, 

k=l 

oo 

n=—oo 

Wormhole conditions for the first order scattering coefficients are: 

(^0 + ^})r=a = e*""(^o + '^l)R=a,c^=e, (20) 
^{^0+^l)\r=a = -e'-^^{^0+^)\R^a,4>=e. (21) 

Similarly for the {k + 1)*'* order coefficients: 
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It is easy to sec that when 1** and A;*'* order scattered waves satisfy wormhole conditions, 
total scattered wave satisfies as well. 



■^\r=a = (*0 + *1 + ^2)|r=a = (^0 + *J + + ^2))\r=a 

k=l 

oo 



fc=l 



Imposing wormhole conditions to (20) and (21), first order scattering coefiicients are 
obtained. 
C-1 yields: 



J2 A„•^(c^a)e^"^ + B^i^W(u;a)•e^"^ = e''*^^( ^^(^«)e''*'+ J2 ClH!^\uja)-e''''), 

%=—oo n=—oo n=— oo 

S^-e-C^ = --|j^(A.-e-^). 
-H A ' [uja) 



C-2 yields: 



K + e^'^^^n = (^n + e-M„). 

-/7«(c.r)|_ 

Solving for 5^ and : 

= -7;r(^«) A + e^'^^Tn (^«)^, (22) 
C;^ = -7+(^a)^ + e---7^(^«)A- (23) 

Note that and are equal to the known parts of (10) and (11), and respec- 
tively. 

A;*'*order scattering coefficients are obtained similarly as: 

= -7„+(c^«)Q' + e'^^^-{coa)Bt (24) 
Cn""' = -7n^(a;a)5^ + e— 7n (^«)C^n- (25) 



IV. 3+1 DIMENSIONS 



In 3+1 dimensions, the solutions of wave equation in spherical coordinates, i.e. spherical 
wave functions, involve spherical Bessel functions and spherical harmonics [13]. In agreement 



12 




FIG. 3: Coordinates used for 3+1 dimensions. 

with the 2+1 dimensional case, ^'o is expressed in terms of usual spherical Bessel functions, 
while and ^2 are expressed in terms of spherical Hankel functions. Referring to fig.3, 

5^ 5^ (26) 

l=—oo m=—l 
00 I 

*i = J^Bim- h\'\ur)Yimie, ip), (27) 

l=—oo k=—l 
00 I 

*2= J] J] a^-/i;'Hu;i?)F;40,<^), (28) 

l=—oo m=—l 

and the wormhole conditions are. 

The addition theorems for the spherical wave functions, for r = c? + ^ are [14], [15]: 

ji{u;r)Yim{9, if) ^Yl c^'rJid)Mu;R)Yi,m'{n - 6, (^), (29) 

I'm' 

h['\u;r)Yim{9, ^) = J2 «/w(^ J^'('^^)^/'m'(vr - 6, <^) for R < d, (30) 

I'm' 

hf\uR)YU'K -Q,ip) = Y,<^\Z'{-d)3v{oJT)Yi,^,{e,ip) for r < d, (31) 

I'm' 

where 
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c(/m|/W|A/i)jA(t^ \x\)Yxi^{x), (32) 



The coefficients c(/m|/'m'|A/i) in terms of 3-j symbols are: 



c(/m|rm'|A/.) = z'+^-^(-ir[47r(2i + l)(2r + l)(2A+l)]V2(^^ ^ J^, (33) 

Tlie expansion (30) and (31) are valid for i? < d and r < d, respectively, and they cover 
region where the wormhole conditions are imposed: R = a and r = a, {a < d). 

Using (29), (30) and (31), ^'c^i and are expressed as a sum of wave functions at 
translated origin as: 



Im Im 

J2 Bim ■ hf\ujr)Yim{e, <p)^Y. Bimjii^r)Yi^{^, ip), 

Im Im 

J2 Clm ■ h\'\ujR)Yim{^, ip)=J2 ClmJl{uJr)Yim{e, if). 
Im Im 

where the analogues of the formulas (5), (6) and (7) are (see appendix A) 

I' 

I' 

Qm = (-i)'+™$]Q'm«i;r(^- (36) 

3-j symbols are zero unless m — m' = fi [16]. Furthermore, d = zd, and Yxf,_{d) = Yxf,_{0, if) 
is nonzero only when = 0. Thus m' = m and that's why the summation over m' drops in 
(34), (35) and (36) 



no(0,^) = t/^, 



o^im'id) = a'^M = E c{l'm\lm\XO)h'^\u;d) '^^ + ^ 



, 47r ' 



c^fm'id) = id) = J2c{l'm\lm\X0)jx{u;d)^ '^^ + ^ 



47r 
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Imposing the wormhole conditions and using the orthogonahty of Yijn{9,(p) for different 
Z,m, yields the 3+1 dimensional analogues of the equations found for 2+1 dimensions: 



givmg 



Bim — e^'^Cim, — 



Im 



Blm + e^'^'^Cim — Q 



ji{uja) 



hf\ua) 



{Clm ~ 6*'^'^ -Sim + -^Im ~ c"^^^im)) 



d_ 

dr 



jiiujr)\r=a 



-{Clm — e^^Bim + Aim — 6*^^^^/^), 



Blrn = -7rI^(^a)C/m + e"^^7„ {uja)Birn + ^im, 

Clm = -'^^{u:a)Bim + e"*'^^7„ (cua)^^ + F;^, 
where Eim and Fim are known functions of Aim '■ 



(37) 
(38) 



Elm = -7n ('^a)Am + e*'^''7n (^o)^;^, 

Fim = -7+(u;a)li„ + e~"^'^7~(u;a)Ai^, 
and '-^^{ujo) are defined similar to 2+1 dimensional case: 




7r(^a) = - I + 



Similar to the 2+1 dimensional case, (37) and (38) can be solved for a ^ 1 and a <^ d 
cases. 



A. Solutions for a <^d and a <^ 1 : 

The asymptotic form of h\^\ujd) for / <^ ud allows us to compute al'^' (d). The similarity 
between 2+1 and 3+1 dimensional cases are remarkable. Indeed for 2+1 dimensional case, 

if we consider Xn = ^k-nHj^\ujd) as an operator on Hn\ojd), the asymptotic form 

fe=— oo 

of Hn\ood) for n <C ud is an eigcnvahic of this operator. Similarly in the passage to the 
3+1 dimensions, considering Xim — (—1)^+"* ^ Xi/m'Oiim ' (d) as an operator on h\^\ujd), 

I'm' 

asymptotic form of h\^\ujd) for I <^ ud is an eigenfunction of Xim- 
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As in the 2+1 dimensional case, the presence of the 7^(a;a) factor at each term of the 
right hand sides of (37) and (38), makes Bim and Ci„i vanish when cua -C /. Thus when 
a <^ d the asymptotic form of h^^\ujd) for / ^ ud can be used. 

For a -C 1, just hke 2+1 dimensions, hl^\Ljd) is zero unless / = and (37) and (38) can 
be solved. 



1. a < d: 

7;'^(a;a) and 7;~(a;a) filter the terms with / > 2uja, thus when a <^ d the only terms that 
contribute to Bim and Q^are / <^ ud. In this case h^^^ {cud) has the asymptotic expression: 



cud 



Then, 



2A + 1 



iujd j 

Bim ^ 5^^Srm'c(rm|HA0)i-(^+^)^ J 

I'm' A 



Substituting 



J I m -m 



c{l'm\lm\XO) = i^'+^-\-ir[An{2l + l)(2r + 1)(2A + 1)]^/' 
gives: 

_ piojd 

Blm ~ J V Bi,y{-ir[{2l + 1)(2/' + if/'dmO 

Ud ^ 

where in the last step the orthogonality property of the 3-j symbols is used [17]: 
Thus, 

Bi^ ~ jV(2^+l) Vi^^'o^V(2^' + l)<^n^o = V(2i + l)r(Ao)(5^o, 

where, for X; being any function of /, the functional T{Xi) is defined as: 

r(x0^5]x,/y(2FTi) 

V 

If m 7^ 0; Bi^ = Ei^, Qm = -F/m and if m = : 

Bio = -e-- — (-l)'v/(2[+l)7r(a;a)r(S,o) + — (-1)' v^(2l + T)7+r(C,o) + £;zo, 
ujcl ujd 

piuid piLod 

Clo = -e-'^^ — {-iy^'{2lTT)-fr{u^a)T{Cio) + — (-l)^vW+l)7/+r(S,o) + F^q. 
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Multiplying each side of these equations by ^'■^/(2^+^) and summing over I gives T{Biq) 
and T{Cio) : 



T(5/o) 
T{Cio) 



1 - e'""^ — TT{{-i)\2l + l)7;-(a;a)) 



AlocL 



-1 



cud 



T(H)^(2/ + 1)7+M) 



■n{-i)\2i + i)7,+M) 



ud 



T((-i)'(2/ + l)7r(u;a)) 



T(i?;o)" 

T(F;o) 



£ a < 1; 

Similar to the 2+1 dimensional case, for a ^ 1, '^f{uoa) becomes a discrete delta function, 
5i. Due to the factors of '-)^{uja) in each term, Bim and are nonzero for only I — m — Q. 
The problem reduces to finding the constants Bqq and Cqo- 



Blm — -Boo (^i 5m, 



I = implies m = and = A, so that 



Soo = ^Soo(5ac(AO|00|AO)/iW(cuo?) 



2A + 1 
47r 



= Booh'i^\ujd), 



Coo = Y,CooSxc{XO\00\XO)h^^\ijd){-iy 



2A + 1 
47r 



-Coo4'Hu;ci). 



i?oo and Cnn are found as: 



^00 



-Boo 




_C'oo_ 





l-e'^^h^Q\ujd) h^^\ujd) 
h'i\ujd) l-e-'^^hi^\u}d) 



-1 






-E'oo 




-^00_ 



B. Multiple scattering 

Multiple scattering formulae for the 3+1 dimensions can be found, by the same steps 
followed as the 2+1 dimensional case The multiple scattering expansion of 3+1 dimensional 
wave functions 



Im 



Im 
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together with the wormhole conditions for the 1** and the A;*'* order scattering coefficients 

lead the formulas for the multiple scattering solution: 

7;+(a;a)Am + e"^''-f^ {u;a)Ai^, 
-f^{u;a)Aim + e''"^''-f~{uja)Aim, 

V. NUMERICAL VERIFICATIONS: 

In this section, the solutions for certain values of a, d, uo and r are evaluated numerically 
for 2+1 dimensions and it is verified that they satisfy wormhole conditions. Numerical 
evaluation of solutions are done by using the multiple scattering results (22), (23), (24) and 
(25). Alternatively (10) and (11) are tested by an iteration method. For iteration, two initial 
test functions B^^ and are picked and substituted to right hand sides of (10) and (11) 
to obtain 5^ and C^. Similarly 5^ and C\ are substituted to (10) and (11) to obtain B^ 
and C^. Continuing this iteration, B'^ and are assumed to converge to the solution. No 
proof for the conditions of convergence is given, it is verified numerically that the solution 
found by iteration method converges to the multiple scattering solution for the parameter 
sets that are considered. 

Moreover, to check the formulas found for a ^ d the solutions found by this method is 
compared with the multiple scattering solution. 

As the velocity of wave is taken as 1 in equation (1), ujd = 27rd/X and ua = 27ra/A where 
A is the wavelength of the wave. Practically if when a light wave in a wormhole universe is 
considered, these values supposed to be much larger (at least order of ~ lO^'^) compared to 
what chosen in the below examples. However, numerical calculations with such large values 
were beyond the capacity of the PC used and there is no reason to think that the formulas 
will fail for large values. 

The incident wave is chosen as a plane wave and A^ — e'"", where a is the angle 
between direction of the incident wave and the y axis. 

Referring to figure 2, the wormhole is located symmetrically with respect to the y axis. 
Consider the reflection operator R with respect to the y axis, i.e. R^{x,y) = '^{—x,y). 
According to the wormhole conditions C-1 and C-2 
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FIG. 4: The contour plot of Rc(^ — e^'^'^i?^) in the vicinity of left wormhole mouth A_. The contour 
circle at r = a shows that — e*'^^i?^)|r=a is constant, {uja = 20; Lod = 120; r = 1; a = tt/3) 




FIG. 5: The contour plot of Re( — + e^'^'^R^)) in the vicinity of left wormhole mouth A_. The 

or 

same contour circle at r = a is evident. (wa = 20; cod = 120; r = 1; a = tt/3) 
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FIG. 6: The contour plot of Re(^') in the vicinity of A_. The incident wave is coming from the left 
with an angle 7r/3 and the shadow is on the opposite side, (wa = 20; cod = 120; r = 1; a = 7r/3). 



-e''^^m>)\r=a = Q (39) 

^^{m + e'-^m)\r=a = ^ (40) 

It is verified that the solution found satisfies (39)and (40) by plotting the contours at the 
vicinity of one of the wormhole mouths. 

In fig. 4, fig. 5, fig. 6 and fig. 7 the parameters are: ua = 20, ud = 120, a = tt/3, ujt = 
1. Figure 4 and Figure 5 show contour plots of the multiple scattering solution for real part 

d 

of "^{x,!/) + e*'^'^i?(\E'(x, y) and —{'^{x,y) — e^'^'^R{'^{x,y)), respectively. In both figures, 

the contour circles at ur = ua = 20 are clearly visible indicating that the values of each 
function are zero along r = a circle. This shows that the wormhole conditions are satisfied. 
The contour plots of imaginary parts -which are not presented here- give the same contour 
circles at r = a. Although the contour is plotted for 0.8a < r < li to make the zero 
contour circle more visible, it should be remembered that the region r < a is not a part 
of the spacetime. Figure 6 is a contour plot of the real part of the solution \I'(r, 0) to give 
an example of a visual image of the solution. Figure 7 is a comparison of the multiple 
scattering solution and the iteration solution. The solid line with markers show the \Bn\ 
that are found by multiple scattering and dashed line with 'x' markers are the difference of 
the absolute values of Bn found by the iteration method and the multiple scattering method. 
The difference is zero for all n; i.e. these two solutions are exactly the same. The results 
are obtained after 20 iterations. The test functions are chosen as constant, = = 1. 
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FIG. 7: Comparison of the multiple scattering and the iteration results. The difference of \Bn\ 
found by these two methods are points with marker 'x' which are zero for all n. (ua = 20; uid = 
120; T = 1; a = 7r/3) 




FIG. 8: Comparison of the multiple scattering and the a <^ d approximation, [ua = 5; cod = 
1600; r = 1; a = 7r/5) 



In fig. 8, the parameters are: ooa = 5, cud = 1600, a = n/5, ujt = 1. This is an example 
ior a <^ d case. In the figure versus n is plotted. The solid line with '+' markers is the 
multiple scattering solution and the dashed line with 'x' markers is the a <^ d approximation 
solutions given by (15) and (16). 

VI. CONCLUDING REMARKS 

The principal purpose of the present work was to investigate the scalar waves in wormhole 
topology. As the wormhole considered is fiat, no complications due to curvature arise. 
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Although the spacetime considered in this work admits CTCs for sufficiently large values 
of time lag r, their existence has no influence on monochromatic waves. The closed timelike 
curves emerge when time lag r is greater than d — 2a. However, r appears in the equations 
only as exp(ici;r). Thus the solution remains the same for all integer /c's where cut — 2kTr + a 
and increasing r does not change the nature of the solutions. This suggests that the presence 
of closed timelike curves does not have a dramatic effect on the scalar wave solutions. 

This should not be surprising considering that, in a wave equation, what really matters 
is presence of closed null curves, rather than closed timehke curves. It is reasonable to think 
that the existence of closed timehke curves will not effect the nature of the solutions as long 
as closed null curves are not present. CTCs are present in the flat wormhole spacetime 
studied here, but still they don't have a significant effect on the solution. The reason is 
explained in [6] : In this kind of spacetimes, the closed null curves are a set of measure zero 
and due to the diverging lens property of the wormhole, the strength of the field is weakened 
by a factor a/2d at each loop in the infinitely looping closed null geodesies. 

The complications related to closed timelike curves are due to difficulty in specifying a 
Cauchy hypersurface when solving the Cauchy problem. Null geodesies are bicharacteristics 
of the wave equation and arbitrary initial data cannot be properly posed in a null direction 
[18]. A spacehke hypersurface never contains vectors in a null direction, thus are good 
candidates for specifying initial data. However there always exist a null direction on a 
timelike point of a hypersurface. In the light of these discussions it can be conjectured that 
no complications arise on the solution of wave equation due to CTCs. The complications 
are mainly due to the nature of Cauchy problem approach. 

To make this point more apparent, let's consider a space-time with topology x S^, 
where is the temporal direction. Clearly this spacetime admits closed timclikc curves. 
In this spacetime, waves from only a discrete set of frequencies can exist due to periodicity 
conditions. But this restriction does not seem to be related to the existence of CTCs, 
because alternatively when a globally hyperbolic spacetime of topology 5"^ x M is considered, 
it is easy to see that it has the same property: Waves only from a discrete set of frequencies 
can exist in this spacetime, cither. This observation also suggests that the restrictions on the 
solutions of wave equation are not due to existence of closed timelike curves. Compactness of 
the topology in either temporal or spatial directions have both similar consequences. Thus 
it is reasonable to expect that waves in a universe where time travel is possible, will have 
similar properties with globally hyperbolic ones. 

If we consider the question in a purely mathematical point of view, the form of wave 
equation considered is almost symmetric with respect to time and space variables. For the 
1+1 dimensions there is complete symmetry (remembering that the minus sign on the time 
derivative does not effect the symmetry since it is always possible to reverse the signs) and for 
higher dimensions the only difference is having more space variables. This suggest that there 
is no strong mathematical background for expecting disparate consequences of existence of 
CTCs compared to existence of closed curves along any space direction. On the other hand 
more space coordinates give rise to asymmetry between the spacelike hypersurfaces and the 
timelike hypersurfaces in Cauchy problem due to the shape of the null cone: Any timelike 
hypersurface passing through a spacetime point intersects the null cone of that point, while 
spacelike hypersurfaces does not. 

There is a strong analogy between 2+1 and 3+1 cases, which suggests that the results can 
be extended to n + 1 dimensions easily. In any dimensions, the solutions can be expressed in 
terms of spherical waves, f{r)Y{fl), where r is the radial distance and fl denotes the angular 
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part [10]. In addition, to be able to apply the same method, an addition theorem similar to 
that of the 2+1 and 3+1 dimensions is needed for this higher dimension. The similarity of 
(10), (11) with (37), (38) suggests that the solution for higher dimensions are readily given 
by these equations where the expressions of B and C in terms of B and C will be found 
using addition theorems of those dimensions. 

Acknowledgement 1 The author thanks Cem Tezer for his valuable guidance and helpful 
discussions, John L. Friedman for his significant comments, and Ozgur Sarioglu. 



APPENDIX A: CALCULATION OF Ai^^ Bi^ AND Cim 

Referring to the equations (26), (27) and (28) 

^Aim ■ ji{ur)Yimi0, <p) ^^Aimji{ujR)Yim{^, ip), 

Im Im 

Bim ■ hl'\u;r)Yi^{9, if) BimJi{u;r)Yim{^, ip), 

Im Im 

YClm ■ h\'\u;R)Yi,n{^,ip) = J2Cimji{uJr)Yi,n{0,<p), 

Im Im 

(29), (30) and (31) can be employed to calculate Ai^, Bi^ and Cim- Considering \E'o, 



*0 = 5^ Aim ■ jliuJr)Yimie, <fi) = Y ^^'n ^ «|w (^IX^^)^'-' (tT - 6, if) 
Im Im I'm' 

= E E c^^jm'{u;R){-iy'^"'%m'{o, 

Im I'm' 

= E(-i)'^'"(E^''-'«l-'^('5)jK^^)i1m(e,(^). 



Im I'm' 

In the last step the order of summations and indices Im and I'm' are interchanged. It is 
assumed that these series converge and changing the order of the summations is valid [15]. 
Thus 

A = {-iy+-{j2Ai'm'a^;::'+{d)). 

I'm' 

Calculation of the Bij^ is identical except 0(\'^''^ is replaced by aj'^' . To find Q^: 
*2 = E ^irn ■ hf\uR)YUQ, ^) = J](-l)'+-Q^ . hf\uR)Yi^{T, - e, ^) 

Im Im 

= ^(-l)'+-am E c^\Z'{-d)ji'{ujr)Yi,m'{e, ^) 

Im I'm' 

= E(E(-i)'^™'^''-'"^-'(-^l)j/K)i^„(^,^) 

Im I'm' 

- E(E(-i)'^'"'^^'-'"^-'(-rfl)jK^r)r,^(^,^). (Ai) 



Im I'm' 



Recalling (32): 



= 5]c(/m|/'m|A0)/.«(a;rf)(-l)y^. 

The ^ Q Q Q ^ factor in c(Zm|Z'm| AO) is zero when Z' + Z + A is odd. Thus (-1)'+' 
and (-1)^ = (-1)'+'', yielding 



Substituting in (Al) 

Im Im I' 

Hence 
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